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ABSTRACT 


There is considerable interest in producing new and improved 
glasses in space by taking advantage of containerless processing. 
However, without gravity-induced movement of gas bubbles, considerable 
problems with bubbles are anticipated. One method suggested for 
fining (bubble removal) is to apply a temperature gradient, so as to 
cause che bubbles to move out because of the dependence of surface 
tension on temperature. Surface tension driven migration of bubbles 
in a thermal gradient has not been demonstrated previously in molten 
glasses. To provide the scientific basis for studying this effect in 
space, ground based experiments were conducted to demonstrate that 
thermal migration actually operated in glass melts. Thermal migration 
consistent with the theory was found in one experiment on a borax melt, 
i.t., there was an approximately linear relation between the bubble 
diameter and bubble velocity for a given temperature and temperature 
gradient. It also appeared that nearby bubbles were attracted to 
one another, which could greatly aid fining. Interpretation of these 
results was not possible because of complications arising from gravity, 
i.e., floating of the bubbles, circulation currents due to buoyancy- 
driven natural connection, and flow of the melt out from the cell. 
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1. INTRODUCTION 

There is one feature common to all aspects of glass melting: 

At one stage or another, raw materials are reacted which either liberate 
gases chemically combined or absorbed in the raw materials and/or entrap 
these gases (p eking voids) in the melting mixture. This leads to bubble 
entrapment. To render the glass transparent and to provide usable 
strength levels, these bubbles are eliminated in a process known as 
'*fining.** While this particular terminology is limited to bulk glasses, 
the requirement for eliminating bubbles is a dominant feature of the 
processing of most glass-containing systems. Ceraniic-to-metal seals 
(high pressure sodium lamps) and glass~to~glass or glass-to~ceramic 
seals (cathode ray tubes) are formed with powdered glass (solder) seals. 
Hermeticity and acceptable strength directly correlate with the elimina- 
tion of the void structure. Predictability and reliability in integrated 
circuits are commonly ensured by the use of passivating glass coatings. 
These are applied by elaborate techniques using powdered glass in liquid 
slurries to assure thin continuous coatings. The typical low expansion 
semiconducting materials require glasses which barely fuse at tempera- 
tures close to the stability limits of some of the components. The 
necessity of eliminating the entrapped bubbles dictates the times 
required at these elevated temperatures - 

Another example of glass processing which is regulated by the 
bubble removal process is conventional porcelain enameling. In this 
process, bubbles must be eliminated to the point of leaving a smooth 
glassy surface. Surprisingly, in this particular application a rather 
well defined bubble structure imparts desirable performance benefits. 

In fact, certain additives, such as clay, are introduced to generate a 
babble structure. The prevention of the defect called '*f ishscale*' in 
sheet iron is one such benefit. 
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1 , I Buoyant Fining 

In spite of the great wealth of experience possessed by glass 

manufacturers, little is understood about the process of fining and 

even less is understood about bubble elimination in sealing glasses or 

1-2 

enamels. The traditional viewpoint of glass fining is based on 

the effect of additives, called **fining agents.” These decompose rapidly, 

to form new bubbles, or the gases diffuse into existing bubbles, 

causing them to grow, whereupon they rise quickly out of the 

melt under the action of buoyant forces. The velocity at which bubbles 

rise in the molten glass may be approximated by the following equation: 

V = (1) 

12n 


where D = bubble diameter 

g = gravitational constant 
p = density of the melt 
n = viscosity of the glass. 

The factor 12 represents a departure from the typical Stokes formula 

3 

because of the nonrigidity of the bubble. 


Simple calculation shows that beneath a particular size bubble, 
buoyant fining does not occur fast enough to eliminate bubbles in the 
times used in conventionaJ processing. For example, a bubble 0.1 mm in 
diameter in a lime glass with a viscosity of 100 p will rise about 
17 cra/day — too slow for commercial processing. An additional mechanism 
is necessary to explain the fining of commercial glasses. 


1 . 2 "Chemical” Fining 

One mechanism for the fining of commercial glasses has been 
elucidated by Greene and Gaffney.*^ Bubbles consisting largely of oxygen 
shrink until they disappear or reach a limiting size associated with the 
fraction of impurity gas originally present in the bubble. Doremus^ 
proposes that the mechanism of fining agents such as As 20 ^ is to decrease 
the concentration of oxygen physically dissolved in the glass and permit 
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greater diffusion of oxygen out of the bubbles and into the glass. 

Carbon dioxide is also able to diffuse out of bubbles into the glass. ^ 

The function of most fining agents can be explained, then, in terras of 
making more oxygen available for diffusion (nitrates), or raising the 
diffusion coefficient of oxygen in the bulk melt (As20^) . These 
mechanisms might be called ’'chemical** fining. 

1.3 "Thermal” Fining 

’’Thermal" fining, to coin a parallel phrase, is a mechanism 
that has been largely overlooked and little understood.^ Temperatur#^ 
differences, with their effect on surface tension, produce a displacive 
force on bubbles, Thermal fining may be the principal mechanism of 
bubble elimination in sealing, soldering, and enamelling operations. 

It is commonly observed that bubbles come to the surface of glass 
enamel quickly regardless of the orientation of the enamel. 

Since seals, solders, and enamels are typically heated externally, it 
is presumed that a driving force toward the external surface is 
established on the bubbles in the fluid glass, 

8“"12 

While significant evidence has been reported which 
confirms that gas bubbles can move in a temperature gradient, no 
evidence has been reported for this phenomenon in glass systems, in 
spite of its obvious importance in glass sealing and enamelling. The 
effect of this mechanism has possibly been observed in experiences with 
direct electric melting of bulk glasses. Contrary to fuel fired furnaces, 
electric melters are hottest at the bottom of the tank. This gradient 
would oppose buoyant fining, and indeed, electrically melted glasses are 
harder to fine, although there are other factors which contribute to 
this . 

While the mechanism of thermal fining is obviously of great 
commercial significance in many high technology applications, it 
represents an even more important mechanism in the space processing of 
glasses. There has been considerable interest in producing glasses in 
a zero gravity environment, particularly because the materials could be 
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produrod without the uso oi a container and its unavoidable contamina- 
tion. In the various equipment schemes proposed to date, no means have 
been proposed for eliminating the bubbles incorporated or exsolved as 
reboil during melting. If artificial gravitational fields were used 
they would nullify the potential benefits of zero G to some of the 
experiments, such as studies of the homogeneous nucleation of 
crystallizing glasses. Migration of bubbles under a temperature gradient 
could indeed serve as a fining mechanism in zero G for the production of 
bubble-free melts. 

There Is a need, therefore, for having a quantitative under- 
standing of the migration of bubbles under a thermal gradient for both 
the control and understanding of ’art-laden* commercial processes and 
for application as a fining mechanism for zero G processing. This 
experiment definition study was therefore concerned with defining 
and developing a scientific basis which would justify experimentation 
in a low G environment and establish the flight hardware requirements . 
Specifically, unambiguous evidence for the migration of gas bubbles 
through a molten glass under the driving force of a thermal gradient 
was sought. The experiments were designed to reduce the influence 
of competitive mechanisms for movement, including surface tension 
driven convection. A prototype hot stage was developed which pointed 
toward a design which would be suitable for flight. The background 
and results of these experiments follow along with an evaluation and 
a proposal for flight experimentation. 
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2, EXPERIMENTAL BACKGROUND 

2.1 Preparation of Candidate Glasses 

The initial set of glasses considered for this experiment is 
listed in Table I, with their composition and batch materials. These 
glasses were mixed in 100 g batches, ball milled 1 hour, and melted in 
alumina crucibles. Each glass was melted several times to determine 
the lowest temperature capable of producing a reasonably well- refined 
glass. Since lead glasses or glasses with high alkali borate contents 
are volatile, a low melting temperature was desirable to remain close 
to the nominal composition and to prevent chemical inhomogeneities due 
to surface volatilization. It was found that glass batches B, C, and 
H melted to good glass with the least difficulty. They were melted 
at 800°C for 30 minutes. Glass I was melted at 900®C for 2 hours. 

Glass I was selected because it was a non-lead glass. These glasses 
were stored in cullet form (chunks) until needed for further experi- 
mentation. 

Further preparation of these glasses consisted of crushing 
and sizing or remelting to introduce bubbles. Once crushed and sized 
the fractions were stored In glass bottles. 

Remelted glass was stirred to introduce bubbles and then 
drawn up quartz capillary tubing or hand-drawn into glass cane, so 
as to produce a bubble array for one of the experimental techniques 
discussed later. 

2.2 Design and Construction of Hot Stage Prototype 

The prototype hot stage is pictured in Fig. 1. The heart of 
the stage is the tapered platinum strip held in position by the two 
copper terminal blocks. Small chrome 1-alumel thermocouples were 
spot welded to the bottom of the strip, allowing the temperature profile 
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TABLE I 

Batch Compositions for Candidate Glasses 


Glass A 

88g PbO 

Glass G 

60g PbO 


21. 4g H^BO^ 


35. 6g HjBO^ 

B 

85g PbO 


17. Ig Na2CO^ 


13. 4g H^BO^ 


lOg AI2O2 


7.5g Si02 

H 

79g PbO 

C 

85g PbO 


19. 6g H^BO^ 


15g Si02 


8g Si02 

D 

70g PbO 


2g AI2O3 


26. 7g H^BO^ 

I 

36g Si02 


15g SIO2 


53g Na2C02 

E 

75g PbO 


24g TIO2 


17. 8g H^BO^ 


12. 4g Li2CO^ 


lOg Al20^ 
5g SiO^ 


7.1g H^BO^ 

F 

70g PbO 
26. 7g H^BO^ 
15g AI2O3 




*The following chemicals were used; 

PbO - Fisher L-71 
” Wisher L-73 

Si02 - Stnd. Cer . Sup. - S925 (-200 mesh ) 

AI2O2 - Fisher A-591 

^®2*"°3 " fisher S-263 

T102 “ Fisher T-315 

L12C02 “ Fisher L-119 
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established by the strip to be determined. Banana plugs connected the 
stage to the power supply and were connected to the copper terminal 
blocks through flexible copper cables. The terminal blocks were mounted 
on a transits base with nylon screws. The sides were brass, £o.<itened to 
the transits base. A simple spring tension mechanism maintained a flat 
configuration of the platinum strip as it heated, thus compensating for 
expansion. The copper terminal block to which it was attached had a 
split base which allowed It to slide freely back and forth in a track, 
while maintaining good electrical contact. 

This stage fit easily on a research microscope sample stage. 

A quartz window fit on the ledge protecting the objective lens from 
high temperature and volatiles emenatlng from the hot glass. A AX 
objective lens with a working distance of 7 mm or more easily focussed 
on the platinum strip surface. 

For the preliminary screening experiments a simple strip 
profile has been used as seen In Fig. 1. This gave an adequate, but 
nonlinear, temperature gradient. R. S. Subramanian of Clarkson College 
has made an analysis of the thermal gradients In the heating strip. 

He has determined the equation which gives the width of the strip, Yj^, 
as a function of axial distance, X, necessary to establish a linear 
temperature gradient (T^-Tj^)/! down the axis of the strip. The equation is 




1^1 


AbL^oT^^ 


1 + BX 

(1 + aX)^ 


( 2 ) 


where 


and 


-'I 

b 

L 

a 


^2 

B 


the current through the strip 

the electrical resistance of platinum at T^^ 

the thickness of the platinum strip 

the length of strip for desired thermal gradient 

Stefan-Boltzmann constant 

temperature at the low temperature end of gradient 
temperature at the low temperature end of gradient 


Pj-Pi 


T -T 
2 1 


with p^, T„ equal to the resistance 
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and temperature of the platinum at the high temperature end of the 
linear thermal gradient. Here Y(X) and X are y(x)/L and x/L, the actual 
dimensions divided by L the length of the linear thermal gradient. 

Figure 2 illustrates how this linear gradient region is incorporated 
into the overall strip. The development of this equation is presented 
in Appendix A» along with the physical conditions : sumed. 




2.3 Preliminary Screening Experiments 

The preliminary screening experiments allowed for familiariza- 
tion with the operation of the hot stage as well as for observing the 
behavior of a molten layer of glass on the platinum strip. It was 
desired to observe the molten layer of glass because there were three 
types of instability that could develop due to horizontal and vertical 
temperature gradients energizing surface tension and buoyant forces. 

As described below, two kinds of cellular convection can be caused by 
a vertical temperature gradient. Whether or not they occur depends on 
the value of the Marangonl number for surface tension and the Raleigh 
number for buoyancy. As indicated by Subramanian's estimate of these 
numbers in Appendix B, they were far below the critical value for the 
molten glass film and, hence, should not have occurred. 

A signii leant surface tension driven circulation was expected 
t^om the horizontal temperature gradient. As derived by Subramanian in 
Appendix C, the maximum rate of fluid flow at the free surface would be: 


U = d/4p dy/dv 
max ’ 


(3) 


where 


d 


film thickness (cm) 


= viscosity (dyne-sec/cm*") 

2 

dy/dx “ horizontal gradient of surface tension (dynes/cm ). If 

the bubble is 1/10 the film thickness, then it is an order of magnitude 
larger and in the opposite direction expected for bubble motiorv which 

= - D. /4p dy/dx 


is estimated to be: 


where D. is the bubble diameter. 

D 


1 


f 


i 



Dwg. 6403A57 


\ _l 

/ 



~r 


10 


Fig. 2 - Schematic of platinum heating strip for hot stage 
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The resultant bubble motion would be very erratic and the bubble motion 
due to a thermal gradient almost impossible to measure. 

Preliminary observation using fine zirconla as a tracer 
showed chat the surface tension flow due to the horizontal gradient was 
indeed active and would probably swamp the effect we wanted to observe. 
This experimental confirmation led us to the modifications described 
in the next section. 

2.4 Experimental Modifications and General Procedure 

The surface tension phenomenon described above was avoided 
by eliminating all free surfaces. Two methods were found to be success- 
ful. The first and most successful was through the use of an inverted 
channel in a quartz window. This method, illustrated in Fig. 3, 
essentially isolated the bulk of the gxass Ircm the small amount of 
free surface remaining. For the inverted channel experiments the glass 
was pulled by hand to make cane 'V'l mm in diameter. A section of cane, 
with Che inverted channel in the quartz window over it, was placed on 
the platinum strip. The strip was heated, melting the glass cane which 
in a few minutes wetted the quartz, flowed into the channel, filled it, 
and established a film of glass everywhere between the platinum and 
Che quartz window. Physical equilibrium was considered to be established 
when the meniscus at the ends of the channel stopped moving. Though 
considerable thermal migration must have gone on during this equilibra- 
tion period, there were still many bubbles present. By raising the power 
to the strip, the motion was intensified (by both increasing the tempera- 
ture gradient and lowering the viscosity). 

The other method used a glass-filled quartz capillary 1 mm 
in diameter. Both eliminated all such surface tension driven flow of 
Che bulk melt. The experimental procedure in each case was as follows. 
The glass to be used was melted, but not fined. If the bubble density 
was not high enough, bubbles were introduced by crude stirring. For 
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Fig. 3- Quartz window with Inverted channel 



the capillary tube experiment, this glass was sucked into quartz 
capillary tubing which was subsequently broken into lengths to fit on 
the hot stage. The experiment consisted of placing the tube on the 
stage and heating enough to melt the entire length of glass. (If the 
entire length of glass was not melted the Internal pressure in the 
capillary caused the bubbles to collapse and disappear.) Because of 
irreproduclbility, it was found to be necessary to improve the thermal 
contact between the capillary and the heating strip by filling the area 
between the strip and the capillary with glass. 

The inverted channel gave the best experimental results on the 
hot stage for several reasons. First, the heat transfer from the strip 
to the glass, which was in intimate contact with it, was much better 
than with a quartz capillary laying on the flat heating strip. Also, 
the small vertical dimension in the channel design allowed almost the 
entire depth to be viewed at one time. This made it much easier to 
follow bubble motion, because the competing buoyancy and surface tension 
effects caused bubbles to also move up and down. With the capillary 
this vertical movement often caused the bubbles to go in and out of 
focus. From the observable hot spot on the strip and the activity of 
the bubbles, it was also apparent that the gradient imposed on the strip 
was much more faithfully impressed on the molten glass with the inverted 
channel than with the quartz capillary. In fact, bubble motion didn't 
occur with the quartz capillary unless the capillary was "stuck" to 
the heating strip with an Intermediate layer of molten glass, and then 
motion occurred as the capillary system "saw" the thermal gradient. 

2.5 Temperature Response and Thermal Equilibration of the Hot Stage 

The thermal response of the hot stage was measured using an 
IRCON Infrared pyrometer. For these measurements, no attempt was made 
to correct for the emlssivlty of the platinum strip, thus all tempera- 
tures were relative. The measured long-term response is shown in 
Fig. 4, while the Immediate response and that due to convection around 
the strip are not indicated. 


Curve 589530 



Fig. 4 -Tomperati'ni equilibration of hot sta»,e.Mea>urernents 
madn with an infrared pyrometer on spots indic'ited 







The first response of the hot stage was to heat quickly, 
within a few seconds after power was supplied to the initial tempera- 
ture shown in Fig. A. After this the temperature drifted slowly up to 
about 50®C. Thi^ drift is thought to have been caused by fir circula- 
tion and/or free convection. Enclosing the platinum strip by putting 
on the glass cover almost completely eliminated this drift. The 
long-term response shown in Fig. 4 was a slight cooling of the "hot 
spot" on the platinum strip and heating of peripheral areas. This much 
slower response reflects the much lower heating rates in the periphery. 

As they heat slowly, more power is dissipated there and less in the 
"hot spot". The result is to cool the "hot spot" and increase the 
temperature of the cooler areas as shown. 

2.6 Typical Temperature Profiles 

As indicated previously, a design equation to enable establish- 
ment of a linear gradient was derived. To cut a platinum strip with 
the precision needed, it would have been necessary to use a flexible 
die. Since the system (and therefore the temperature and the gradient 
that was best to use) was not yet established, a carefully profiled 
heating strip was not used in these screening experiments. The 
alternative used was a simple tapered design like that indicated in 
Fig. 1. This produced a hot spot and a steep temperature gradient in 
both directions toward the copper posts. Thermocouples were spot 
welded to the underside of this platinum strip and from these the tem- 
perature profiles shown in Fig. 5 were determined. The highest tempera- 
ture profile in Fig. 5 is the type used in our bubble migration experi- 
ments. Unfortunately the strips with thermocouples welded to them were 
too badly warped to be used successfully for the actual experiments. 

It is from this measured profile and others that we estimated (conserva- 
tively) the 300®C/cm gradient that is used for the calculations in 
Section 3. 
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Thermocjuple Locations 



Thermocouples Positioned at .32 cm Intervals 


Fig. 5 - Measured temperature profiled for 
platinum strip shown in figure 1 
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3. RESULTS AND DISCUSSION 
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3.1 Thermal Migration in Common Lead Glasses 

It is generally assumed that the surface tension of a liquid 

13 

always decreases with increasing temperature. With this assumption 

in mind, the glasses B, C, H, and J, were used for the initial bubble 

migration experiments. The negative (migration toward lower temperatures) 

and/or contradictory results obtained using these glasses was puzzling. 

Surface contamination or some other experimental problem was suspected. 

An examination of the literature brought to light the work of Shartsis, 

14 

et al who measured surface tension in several binary silicate and 

borate systems. The estimated surface tension for glasses B, C, and H, 

based on their data, is given in Table II. The three glasses each show 

14 

Increasing surface tension with increasing temperature. Though not 
large, this explains the contradictory results observed in the experi- 
ments using these glasses. This did not eliminate the possibility of 
other effects, but did indicate that the originally expected motion 
toward higher temperatures should not be observed. 

The surface tension temperature coefficient of glass C in 
Table II is positive and about one-half the magnitude of borax (see 
next section). This still should have produced a readily observable 
buLLle motion, but motion was only observed sporadically. This may be 
explained as follows. With a positive variation of tension with tempera- 
ture, both thermal migration and buoyancy would act in the same direction 
causing bubbles to move upwards. This is what was observed. Bubbles 
rose, coalesced and were pinned at the surface, showing only halting 
movement at best. The combined buoyant and thermal gradient forces 
pinned bubbles to the top of the channel where any surface irregularities 
would interfere with lateral motion. Bubbles that did not finally work 
their way into the high gradient region drifted outward, but the motion 
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TABLE II 

Calculated'*’ Surface Tension for Lead Silicate and Boroslllcate 
Glasses. Based on the Experimental Data of Shartsls, L., 
Spinner, S., and Smock, A. W. , J. Research Natl. Bureau of 
Standards, 61 (1948). 





Surface Tension (dynes /cm) 


Glass (wt%) 


600“C 

700*C 800“C 900* C 

1000 "( 

B - 85 - 

PbO 




7.5 - 

®2°3 

- 

171.8 174.7 177.5 

- 

7.5 - 

S 102 




C - 85 - 

PbO 


182.5 187.3 192.4 

194.6 

15 - 

S 102 




H - 79 - 

PbO 




11 - 

®2°3 

_ 

189.0 189.7 189.9 

— 

8 - 

S 102 




2 - 

AI 2 O 3 * 




J - 73 - 

PbO 




11 - 

®2°3 

_ 

_ « 


11 - 

AI 2 O 3 * 




5 - 

S 102 





“^Additivity of surface tension for the PbO>S 102 and PbO-B^O^ systems 
Is assumed. 

Effect of alumina In these systems Is not known, therefore, no attempt 
was made to calculate values for glass J. 
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showed the effects of interference from the roof of the channel. This 
system could only be studied effectively in a low G environment. 

3.2 Thermal Migration of Bubbles in Borax Glass 

Thermal migration of bubbles in molten borax was observed in 
both a quartz capillary tube and with the inverted channel technique. 

The best observations were made with the inverted channel because its 
shallow depth (500 Mm) allowed all bubbles to be seen at once and 
because there was good thermal contact with the heating strip, as 
discussed previously. 

The thermal migration experiments were performed using the 
prototype hot stage. It was placed on the stage of a Zeiss Ultraphot 
III-B microscope* The heating strip was illuminated with the microscope 
lighting system and viewed in reflected light using a 4X objective lens, 
giving a total magnification of SOX. The working distance of the 
objective lens allowed the quartz cover of the hot stage to remain in 
place, as required for good temperature stability of the heating strip. 

The inverted channel was filled using another technique. 
Bubble-containing glass rods were hand drawn from a crucible melt of 
borax. Rods of about 1 ram diameter were produced. A section of this 
rod '^2.5 cm long was set on the cold heating strip. The quartz window 
containing the inverted channel was placed on top of the rod so that 
the rod was in the channel. The system was heated to completely melt 
the borax glass which then flowed wetting the quartz and filling the 
channel and the narrow seam between the rest of the quartz window and 
the platinum heating strip. The final equilibrium condition is illustrated 
in Fig. 3. 

The bubbles, particularly in the hottest section, moved 
and coalesced during (mechanical) physical equilibration of the system. 

This resulted in a f«*w larger bubbles in the center of the .lOt zone; but 


19 


in general, most of the bubbles in the cooler section did not have time 
to move very far in the first 1 to 2 minutes. The system was considered 
to be in mechanical equilibrium when bubbles trapped between the platinum 
heating strip and the quartz window (not In the channel) stopped moving. 

At this point several photographic sequences were taken with time Intervals 
of 5 or 10 seconds in the sections with rapidly moving bubbles, and 20 
seconds In the outer region with slower moving bubbles. It was noted 
that the bubbles trapped between the platinum and the quartz window 
(outside of the channel) showed little or no movement. Bubbles In the 
channel generally moved up the thermal gradient as was expected. 

The bubble migration observed is well illustrated by Fig. 6. 

Here bubbles next to the channel wall (dark band just above two larger 
bubbles) are shown migrating toward a region of higher temperature. 

Two features to look for in this series of pictures are: 1) the 
acceleration as the bubbles move to higher temperature, and 2) larger 
bubbles move faster than smaller ones. The curved line defined by the 
bubbles in the series of eight pictures is due to acceleration. The 
difference in velocity due to size is Illustrated by the second largest 
bubble Overtaking the smaller one just ahead (to the right) of it, by 
the 20-second frame. Of significant Importance Is an apparent attraction 
of bubbles to one another, enhancing coalescence. This phenomena has not 
been reported before. 

3.3 Comparison of Experiment with Existing Theory 

The motion of an isolated bubble in a liquid with constant 
properties is given by the approximate formula of Young, Goldstein and 
Block (1959); 

*-D/4n dy/dx (4) 

where v^ * velocity of the bubble 
D ® diameter of the bubble 
n * viscosity of the liquid (melt) , and 
dy/dx ® surface tension gradient. 

(See Appendix D.) 
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Fig. 6 - Bubble migration in a thermal gradient in molten 
borax. Large bubble has a diameter of 260 urn. 
Note bubbles accelerating as they move to higher 
temperatures. 
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In Fig. 7 a series of curves were plotted which correspond to 
the dependence of on the bubble diameter, D,for a sequence of tempera- 
tures. Values of the viscosity, n, were Interpolated from the reported 
values of Shartsls et al^^ for the Na20-B20^ system. Their data (two 
points) were fit with the curve log n ■ 10.221 - 12153/T(^C). The 
surface tension gradient, dy/d\, was calculated from the data of Shartsl3 
and Capps ° giving dy/dT - -.076 dynes/cm-°C and a reasonable estimate of 
dT/dx of 300**C/cm, giving dy/d\ » -22.8 dynes/cm . The data points 
plotted In the figure are bubble velocities determined from photograph 
sequences such as the one shown In Fig. 6. Three sequences were taken 
at approximately the same gradient and temperature conditions In different 
experiments. The data agree with the prediction that bubble velocity 
Is proportional to bubble diameter. The dashed line shown in Fig. 7 Is 
a visual fit to the observed bubble velocity-diameter trend. Using the 
velocity equation, a temperature of 785**C Is estimated to have been the 
temperature at which these bubble velocities were measured. Data points 
for bubble velocities In a cooler section of the channel Ind.tcate a 
temperature at or below 700°C, which again Is consistent with the tempera- 
ture estimated for that region. 

There are several conditions that might account for the scatter 
of the data points: 

* There was a steep vertical gradient present as well as a 
horizontal gradient. Therefore small bubbles on the bottom saw the 
highest average temperature, and would have moved most rapidly. 

* The vertical size of the channel was '^'500 ym. Large bubbles 
(<250 uq) , therefore, were always within one radius of a boundary where 
any Interaction would have slowed them down. 

* There was some Indication that bubbles tended to stick to the 

walls. 
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Bubble Diameter (|im) 


Rg. 7- Comparison of predicted and measured bubble movement 
rates in molten borax. 

o " measurements at high T near center of Pt strip 
measurements at low T 

Solid lines calculated from approximate equation 
assuming dT/dX = 300®C/cm. Dashed line drawn to fit high T 
data, with T ~ 785, ®C estimated by interpolation 












* Bubbles apparently attracted each other, at least part of 
the time, so that bubbles separated by a distance less than about a 
diameter tended to interfere with each other. This was noted in measuring 
bubble velocities. Large ones appeared to drag smaller bubbles or even 
cause them to reverse their motion. 

* Buoyant forces acted on the bubbles. 

The fact that all the bubbles were not at the top of the 
channel in these experiments was probably due to the strong vertical 

g 

gradient driving .:he bubbles down as in the experiments of Young, et al. 

If dy/dt were not as large as it is for borax, or of opposite sign, the 
bubbles would all have been at the top of the channel. 

3.4 Approximate Comparison of Thermal Gradient 
and Buoyant Forces in Molten Glass 

The velocity of bubbles rising in a molten glass due to buoyancy 
is approximately |v| = ♦ Deviation from Stokes formula is due to 

nonrigidity of the bubbles. The velocity of a bubble due to a thermal 
gradient given previously is |v| = D/4n dy/dx. Equating these velocities 
(equivalent to equating the forces) gives 

P ^ 3 dy/dT dT/dy (5) 

(cm) gp 


- .0001 dT/dx (®C/cm) for borax 

where dy/dT = .076 dynes/cm-^C 
p = 2.4 grams/cm^. 

In other words a bubble 300 pm in diameter in molten borax would not rise 
in a vertical doimward thermal gradient of 300*’C/cm. Smaller bubbles 
would move downward. This type of phenomena was observed In the experi- 
ments using moltpu borax. 
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4. SIGNIFICANCE OF RESULTS 

For the first time in molten glass, the migration of 
bubbles driven by a thermal gradient has been obseirved. This mechanism, 
while perhaps commonly known in other systems, has not been identified 
as fundamental to glass processing technology in any of the current and 
classic texts. The results of this work show that for bubbles 300 ym 
and smaller and for thermal gradients which can reasonably be expected 
in sealing, enamelling and glass encapsulation (100-300*C/cm) that 
thermal fining dominates buoyant fining. 

Even at this preliminary stage, these observations have 
implications vrtilch bear on the techniques of sealing, enamelling, 
encapsulation and passivation. They point to the need to generate the 
largest thermal gradients through the bubble-containing glass thicknesses 
as can be tolerated by Che rigid components. The long-standing approach 
of substituting tlme-at- temperature for peak temperature excursion can 
only be beneficial Co specimens oriented so as to be beneficted by 
buoyant fining. 

Composition is seen to be a critical factor in the magnitude 
of the surface tension force. It has been seen that in the commercially 
very important Pb 0 -B 20 j-Si 02 system, a striking dependency of the surface 
tension temperature coefficient is Co be found. Slight alterations in Che 
PbO-B^Oj ratio swing this coefficient from negative to positive through 
zero.^^ The substitution of Si 02 for the ^ 2^3 to drive the coeffi- 

cient positive. These effects can have significant implication on the 
design of passivating glasses. Undoubtedly in the past^ trial-and-error 
has led Co Che effective optimization of compositions to maximize thermal 
fining in many commercial sealing and enamelling systems, without a conscious 


25 


knowledge of this effect. With a firmer understanding of composition's 
role on surface tension coefficient, better designs could be produced 
more quickly. 

The evidence of these experiments provides a high level of 
confidence that bubble migration in thermal gradients (thermal fining) 
is sufficient to serve as a method for the elimination of bubbles in 
glasses prepared in space. As pointed out, with reasonable gradients, 
bubbles could be driven out at rates comparable to those achieved on 
earth with buoyant forces. In order to apply this technique in the 
space environment, a much better understanding of the process must be 
obtained. This has been impossible under gravitational constraints. 

The difficulties are elaborated in Section 5. 
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5. NEED AND JUSTIFICATION FOR LOW-G EXPERIMENTS 

We have demonstrated that a temperature gradient can cause 
bubbles to move in a temperature gradient, due to the dependence of 
surface tension on ten^)erature. However, for several reasons It Is 
lxq>osslble to do glass fining experiments on earth which are free of 
the effects of gravity. Most obvious Is buoyancy-driven natural con- 
vection, which can cause rapid global movement of the glass melt and 
Irregular thermal losses from the apparatus. Free convection In the 
melt can be avoided by use of a vertical geometry, taking care to avoid 
horizontal temperature gradients. Convective heat losses In the 
surroundings could be avoided by operating In a vacuum, for example. 

The other gravity Induced problems outlined below can only be avoided 
by performing the experiments In space. 

If a single bubble Is used In a vertical apparatus, then 
bubble movement rates can be measured. However, because of buoyancy 
of the bubble, the movement rate will be quite different than In space. 
Nevertheless, such ground-based experiments can be used to test approxi- 
mate theoretical results. For example this has been done by Young et al. 
to test their theoretical result. In which It was assumed that the 
movement rate was zero, l.e«, the temperature gradient was adjusted so 
as to exactly compensate for buoyancy. 

In glass fining we have the added complication that many 
bubbles are present, with a variety of sizes. Thus with gravity they 
will all rise at different rates, with the large ones overtaking the 
small ones. In addition, as we have observed experimentally, the 
bubbles seem to be attracted to one another In a temperature gradient. 
Such behavior can only be studied unequlvocably In space. 
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5.1 Theoretical Work Needed in Support of Rocket Results 

g 

As noted earlier, the theory of Young et al incorporates 
several assumptions: 

• A single bubble in a large volume of liquid. 

o Very low movement rates (actually Vg -> 0) . 

• Constant properties . 

o No volatile components. 

None of these assumptions are valid for glus^ ining. Merely doing 
experiments is insufficient to establish quantitative methods for proper 
design of ^pace processing of molten glasses. New theoretical work is 
absolutely essential. In order of increasing difficulty this theoretical 
work should treat: 

1. Finite movement rates without gravity. 

2. Influence of variable properties, especially dependence of 
viscosity on temperature. 

3. Influence of volatile component. 

4. Interaction of bubbles with one another and with bounding 
surfaces of melt. 
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6. PROPOSED ROCKET EXPERIMENTS 

Three rocket flights are proposed. The purpose of the 
flights are to obtain data in the absence of gravitational interference 
which will permit refinement and confirmation of the theoretical model 
for thermal fining in glasses. These experiments must encompass varia- 
tions in surface tension, surface tension temperature coefficient, 
temperature, temperature gradient and volatility to develop an adequate 
theoretical model. 

6.1 Glass Systems 

Two glass compositions are proposed for evaluation. The first 
is the Borax glass, sodium tetraborate, and the second is a PbO-SiO^ 
glass identified as glass C in Table I. The first glass is probably 
typical of most glasses in that it has a negative surface tension 
temperature coefficient. Two experiments are proposed with this glass. 
Bubble array motion would be studied with temperature gradients targeted 
at 100 and 300®C/cm, at an average temperature of approximately 780®C. 

The bubble arrays afford the advantage of being able to record data 
from a large number of bubbles of different sizes. This has the obvious 
advantage of improving the statistics of the data collection. In 
addition, apparent bubble interactions can be observed and these effects 
analyzed. 

The third experiment will study the motion of bubble arrays in 
the PbO-SiO^ glass in an intermediate temperature gradient (approx. 200®C) . 
This system has a positive surface tension temperature coefficient causing 
bubbles to move to the colder end of the melt. 

Where the precision of new theory dictates, measurements of 
surface tension and viscosity versus temperature will have to be made. 


29 








¥ 


! 


I* 


1 



6.2 Data Treatment 
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6.2.1 Theoretical Model 

Theoretical work will be carried out by Clarkson College 

g 

personnel to derive an improved model over Young, et al which will 
attempt to treat variable viscosity temperature properties, bubble 
interactions and component volatility. 

6.2.2 Data Analysis 

All of the photographic data will be converted to double notched 
16 mm film and analyzed on the Clarkson College motion analyzer, by 
Clarkson personnel. This will permit accurate measurement of bubble 
position on each frame and automatic recording on computer cards. 
Telemetered temperature data will be converted to temperature gradients 
vs time. Results will be compared with the new theories. 

6.3 Flight Hardware 
6.3.1 Hot Stage 

The hot stage will be one of two designs. Primary considera>- 
tion will be given to reinforcing the design of the successful strip 
heater which uses the inverted channel and quartz cover slide (Fig. 1), 

A thicker platinum strip will be used to minimize the deformation around 
the thermocouple connections. In order to reduce temperature gradients 
through the thickness of the channel the back surface of the cover glass 
will be coated with an IR reflecting film of doped tin oxide. This is 
sufficiently transparent at these temperatures not to interfere with 
the optical measurements, yet it is thermally stable and an efficient 
IR reflector. The shape of the strip will be designed as in the past, 
but actual temperature profiles will be measured with welded thermo- 
couples. 
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An alternative design will be briefly and initially considered. 
This design will simply consist of a small diameter fused quartz tubes 
with two Independent heating colls attached to each end (Fig. 8). A doped 
tin oxide film would also be used on the outer surface of this tube in 
the gradient region to reduce transverse temperature gradients. The 
operation of this arrangement would be somewhat different in that both 
ends would first be heated simultaneously to the same temperature (the 
cold temperature) , and then one end would be heated to the hot tempera- 
ture to establish the gradient. The one thing that needs to be established 
in a simple preliminary earth-based experiment is to see that the equili- 
brium gradient can be established before all bubbles have migrated to 
one end. If this is the case, this simple design will replace the strip 
heater. Flight times may require treating transient gradients. Since 
multiple tubes can be used, data on several systems can be obtained in 
each flight. 


6.3.2 Camera Requirement 

Ideally photographs would be taken at a rate of one frame/s 
or greater. Continuous tungsten illumination has proven adequate. 

6.3.3 Microscope 

Simple, low power reflected light microscopy has proven to be 
sufficient in the ground based experiments. Either the reflectivity of 
the strip heater or the bubbles themselves has produced adequate contrast 
for good photographic coverage. Enhancement techniques such as polarized 
light and Nomarski phase contrast techniques did not prove necessary. 
Magnification of 20 to 40X would be ideal. 

6.3.4 Controls 

Necessary controls and power conditioning for the illumination 
and hot stage will have to be designed and constructed. Choice of 
techniques will be dictated ly vehicle requirements. 
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Dwg. 6413A64 



Fig. 8 ~ Alternate cell design for bubble migration 
experiment. Linear gradient through the cell achieved 
by controlling the temperature in the end caps. 
Thermocouples at each doped tin end of the cell used to 
control power to end caps. Oxide film reduces heat 
loss through the surface of the cell while allowing 
good optical observation 
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6.3.5 Testing 


I 


The system components will be tested to assure proper operation 
individually and as an integrated system under normal ground based condi- 
tions. Although the system will be designed to meet the anticipated 
flight environment, flight qualification testing will have to be negotiated. 

6.3.6 Integration 

The system will require a NASA-furnished multiple exposure 
camera and electrical power at 28 V (nom.) dc. Power requirements would 
fall in the range of 100 to 500 W, depending on the choice of illumination 
and hot stage control design. Camera interface requirements, power/heat 
dissipation constraints, and envelope limitations are required from NASA 
in the early phases of the equipment design programs. 

6.4 Expected Results 

6.4.1 Short Term 

As a result of the proposed experiments, new and reliable 
theoretical models should be obtained for bubble migration due to 
thermal gradients in a molten glass. These formalisms should prove to be 
of great and wide-spread interest to the community of glass technologists. 
The implication of these findings should be felt in the design of glass 
compositions for semiconductor passivation, encapsulation, glass solder 
sealing, enamelling, and direct electric melting. 

6.4.2 Long Term 

The theoretical treatment of bubble migration in glasses due 
to thermal gradients would serve as the foundation for the adaptation 
of this mechanism to the fining of glasses, ceramics and metals manu- 
factured in space. Melting techniques, post heating techniques and 
equipment will be designed, based on these results, which will be able 
to generate the indicated thermal gradients to fine space-melted materials. 
This will prove to be a powerful adjunct technique to the basic complement 
of space processing equipment. 
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APPENDIX A 




Design of the Heating Strip for Glass Fining Experiments 

The purpose of this exei^cise is to design the shape of the 
Pt-lO/Rh heating strip so that a desired temperature gradient may be 
produced along its length* In keeping with the spirit of simplicity 
in these preliminary calculations, several assumptions will be made — 
these may be justified on intuitive grounds and also supported by rough 
numerical calculations* 

1* Steady state conditions will be assumed* 

2* Under these conditions, the principal heat transfer out 
of the strip occurs by radiation from the bottom surface of the heating 
strip. There will be a certain amount of heat lost through the molten 
glass from the top and sides* However, an attempt to take this into 
account will lead to modeling which would be too complex to undertake 
at this stage and may not be worth the return* (Probably, this will be 
weakest assumption.) There will also be heat lost by radiation from 
the sides of the platinum strip which, for similar reasons, will be ignored. 

3* The bottom surface of the heating strip will be assumed to 
be a perfect black body* This may be achieved approximately by coating 
it with platinuia black* 

4* Temperature variations over the thickness of the heating 
strip will be Ignored — it can be shown that such variations should not 
be more than a fraction of a degree Celsius. 
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It may be seen from the physics of the problem that the heating 

strip should be symmetric about the X-axis. The design question is 

9T 

"What is yj^(x) so that “ 8 where g is a specified constant?" Under 
the assumption of negligible temperature variations in the z-dlrection 
(normal to the plane of the paper), the differential energy balance for 
the heating strip takes the following form: 


kVT - i cI* + - 0 


( 1 ) 


The symbols are defined below: 


V 

T 

b 

a 


* Two-dimensional gradient operator 
■ Temperature in the strip in Kelvin 

T =• T (x,y) 

» Thickness of strip, cm 

* Stef an- Boltzmann Constant 


5.67 X 10 


-12 


W 

ci?i?» 


P 

k 


« Current through the strip, amps (it has been assumed here 
that the current distribution in the y-directlon is uniform; 
this is not strictly true, but is a good approximation for 
small yj^(x). The level of approximation introduced by this 
assumption is probably consistent with that of the rest of 
the calculations* as will be seen later). 

« Specific resistivity of material of strip; ohm-cm,p « p(T) 
s* Thermal conductivity of material of strip; 


For the physical problem under consideration, one may write the following 
set of boundary conditions. 

T (o,y) - (2) 

(3) 


T (L.y) = T 2 
0 

7T • n ■ 0 on y 

Where "n is the normal to the curve y 
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(4) 

(5) 


y^Cx). 
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Equation (.5) may be rewritten as 
9T , , V ^ 3T ^ 

■ ^1 37 ‘ ° y 


y^Cx) 


( 6 ) 


Normally, y^|^(x) is specified and the problem for T(x,y) is completely posed 
and may, in principle, be solved. However, here we wish to determine y^(x) 
given the requirement that 

T(x,y) “ + H (y)x (7) 

It is quite possible that no solution of the form of equation (7) exists 
for equations (1) thru (6). Anyway, let us Introduce it into equation (1) 
and examine the result. Since it is known that the thermal conductivity 
of Pt - 10% Rh is practically constant over the temperature range of Interest, 
the result is 


kxH" - a (T^ + Hx)^ + 




(8) 


It is clear from equation (8) that a solution of the form of equation (7) 
is not suitable unless H is a very weak function of y. In this case, 
equation (8) may be satisfied by the approximations 


i (T^ + Hx)^ 2 


4b-^ 


H" * 0 


(9a) 


(9b) 


where H is a transverse average of H(y). In effect, equation (9a) strikes 
a balance between elec ri cal dissipation and radiation from the surface 
ignoring conduction in the x-direction (which is being taken into account 
because of the change in cross-sectional area) . Does equation (9b) have 
a meaningful solution? Since H'(o) = 0, the only solution is 

H = constant (10) 

The boundary condition at y = y^^(x) gives 

-Hi y,'(x) + H'| x = 0 (11) 

iy.y^ 'yr/i 

with H * constant, this can be satisfied only if yj^(x)~ o. 

What we have learned from the above is that a linear axial temperature 
gradient can be achieved only approximately. The approximation would 
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become better for smaller y^. So» our best hope Is to use equation (9a) 
(which could have been obtained from a simple one-dlmenslonal energy 
balance assuming no dependence of T on y) to obtain the required shape 
y^(x) and examine y 2 ^(x). The result is 


yI (x) 


4bof 


iffi. 

/tT** • 


(12) 


For a linear tenq>erature gradient, we may write 

T - Tj^ + (T 2 - T^) ^ (13) 

or e « T/Tj^ = 1 + a X (14) 

where X “ T ^ . If it is assumed that the specific 

L Ij 

resistivity p is linear in temperature over the range of teiqieratures 
involved. 


or 


P •’ PI + 

♦ = P/Pj 


P2-P1 

t7It7 ^’^“’’1^ 


= 1 + ex 


p + (p “p ) T 


(15) 

(16) 


where 

Defining 


g = , [ Here p 

Pi 


Y “ ^ so that Yj^(X) 


j=p(Tj^) and p 2 =p(T 2 )] 


yJ (X) 




. JL±Jx 

^ (l+aX)** 


} 


Since y^(x) = Y^(X), all we need to do is to perform a sample calculation 
with representative parameters to judge the degree of approximation made 
here. Tu; following values are chosen as typical for the parameters. 


T, = (450"C) - 723*K 

^ Pq®c'’ 

z (550“C) = 823"K ^ 

p(T^) - pj - 1.725 p(t - 0-C) 


22.45 pn 

=- 3.87 X lO"^ £l-cm 
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p(T,) - p » 1.88 p(t =■ O'C) = 4.22 x lO”^ n-cm 


0.138 


8 “ P - i z 
Pi 


C4.22 - 3.87) X 10“ 
3.87 X 10“^ 


- 0.09 


b 2 2 mils = 5 X 10 


we here 


o “ 5.67 X 10 

1=8 amps 


Y^(X) 


-12 W 
cm^K^ 


, 1/2 (1 + 


and» from this. 


Yj'(X) 


(1 + aX)' 


e-4a-3agX 

2(l+aX)^(l+eX)^'^^ 


Here, the dimensionless group A is given by 

A * ' 4bL^ I ^ * 0e08 for the above values of the parameters* 

The table below gives values of and Y^*(x) calculated at various X-values 
using the above results. 


X 

Yi(X) 

Yi'(X) 

0 

0.283 

-0.0653 

0.1 

0.276 

-0.0629 

0.2 

0.270 

-0.0606 

0.3 

0.264 

-0.0585 

0.4 

0.259 

-0.0564 

0.5 

0.253 

-0.0544 

0.6 

0.248 

-0.0526 

0.7 

0.243 

-0.0508 

0.8 

0.238 

-0.0491 

0.9 

0.233 

-0.0474 

1.0 

0.228 

-0.0459 
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The shape of the strip is plotted in Figure Al. 

Even though the shape of the strip appears almost trapezoidal. 

It is important to adhere to the exact shape given by equation (18) instead 
of a rough trapezoidal design since the table reveals a 42Z change in 
(based on the smallest value). From the physics of the problem it may be 
seen that the strongest variation of temperature with the y-coordinate 
would occur near the edges y =±yj^(x). Therefore, it is important to focus 
the microscope in a region in the center of the strip around the x-axis. 
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Fig. A1 - Sample design for Pt-10% Rh heating strip 








ADDITIONAL ANALYSES OF TEMPERATUKL VARIATION 
WITHIN THE STRIP 

Temperature Variation Across the Thickness of the Strip 

The maximum rate of heat loss from the strip due to radiation 
may be calculated as 

ot/ :: 2.6 ^ 

An upper estimate of the temperature gradient in the z>direction may be 

obtained by assuming that all this energy is generated at the top surface 

of the strip and conducted through its thickness to the radiative surface. 

W 

The normal conductivity of Pt-10% Rh is ~ 300 “K (and is probably 

larger at a higher ten^.) which gives a 3T = 2 .6W/cm^ oov/ 

3z 0.3W/cm‘K ^ 

-3 

Since the thickness of the strip is 5 x 10 cm, the temperature variation 
across the strip is AT s 0«05^K and may be safely Ignored in this analysis. 


Temperature Variation in the y-Direction 

It appears that the temperature will decrease with distance 
from the x*axis in either direction so that the effect of temperature 
variation with the y-coordinate will be to force the bubbles to migrate 
towards the centerline. The maximum I—) will occur at the boundaries 
y « ±yj^(x) and a rough estimate may be obtained from the boundary condition, 
equation (6)« 

--Si 


Max 


3y 

il 

3T 

3x 


Max I yj'(x) | -Max |y^'(X)| - 0.065 


The actual ratio is probably smaller. In other words, the transverse 
temperature gradient will be much smaller than the axial temperature gradient. 
Such a gradient can probably be tolerated. It is iiq>ortant that the heat loss 
from the edges y - ±y^(x) be held to a minimum since the effect of such loss 
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3T 

would be to Increase |~| and promote transverse variations in temperature. 
One way to acconq>lish this is to polish the edges to reduce radiation 
loss (which will be the principal means of heat loss in SPAR experiments). 


Validity of the Assunqption of Negligible Heat Loss From 
The Strip to the Glass Film 

At steady state, with proper radiation shielding of the quartz 

cover on the glass film, there should be negligible radiation loss from 

the top surface of the quartz cover. So, the maximum heat loss from the 

strip to the glass film and quartz will probably be on the order of 

magnitude of the conduction rate in the glass film (and quartz) in the 

3T 

x-directlon. Assuming an approximate gradient 100 average 

thermal conductivity of 0.02 — ^ for the glass film and the quartz 

cm iC 2 

cover » and a total cross sectional area for conduction < 0.05 cm (from 
glass film depth == quartz thickness » 0.025 cm and width s 1 cm)^ the 
heat conduction rate in the x«direction is roughly 


W 


q ^ 0.02 


ov X 0.05 cm X 100 — 
cm- K 


% 

cm 


= 0.1 W 


The heat lost from the bottom of the strip through radiation is 


q > 1.3 W 


assuming a minimum width of 1/2 cm for the strip. Thus^ the heat loss by 
radiation will be the dominant factor balancing electrical dissipation in 
the strip. 
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APPENDIX B 1^’ 


Surface-Tension Driven Flow in the Molten 

m a ee V-dm ui ^h a Prtf»A Sii'rfilfA 



We shall make several reasonable assumptions in order to get a 


first estimate of the nature of the flow to be expected. 

(1) Flow is assumed to be steady — the fluid is bounded by 
suitable barriers at x==0 and L. 

(la) The length and width of the film are much larger than the 
depth. Therefore, variation in the z-direction can be ignored (3/3z 5 0), 
and away from x*0 and x*L, 3u/9x 5^ 0. Here, u is the x-component of the 

‘ velocity field. 

(2) Density variations with temperature will be ignored. 

I Buoyancy driven convection may be ignored on earth as a first approxima- 

I 

^ tion. This is reasonable in view of the large surface to volume ratio 

I in the film. 

I (3) Even though the surface will be deformed (for a balance 

? 

I of normal stresses), it will be assumed flat for convenience. The 

I 

i equation of continuity gives 





f I 


Here v is the y-component of the velocity. Since v(0)»0, v*0 everywhere 
in our region of interest* Therefore, the variation of pressure in the 
y-direction is essentially hydrostatic and may be ignored for convenience. 
The principal variation of the pressure p will be in the x-direction 
under the assumptions, the only equation of motion that we need is the 
one for u. It simplifies to 

i- .. in = n\ 


The boundary conditions on u(y) are: 
no slip at 

the wall “<“> ■ 


balance of 
tangential stresses 
at the surface 


(Here 'JT is the surface tension.) 


While |i will depend on y to a certain extent because of transverse 
temperature variations, its principal variation will be with x. So, it 
may be treated as a constant in the y-integration. The resulting 
velocity profile may be obtained by a straightforward integration of 
equation (2) to (4) to give 


2\i dx 


2 yd) + ^ ^ 
y dx 


From the condition of no net flow at any axial location. 


VA 

j u dy 


the pressure gradient may be related to the surface tension f radient. 
Using equation (5) in (6), the final result is 


t 


I 


(7) 


I 


dx 


1 _ ^ 

2d dx 


This result may be used in equation (5) to rewrite the velocity field 
completely in terms of the surface tension gradient. 


u(Y) = f- ^ - 2Y) 

4 VI dx 


( 8 ) 


Here Y = y/d. The dimensionless velocity profile, U(Y) versus Y, is 
plotted in Figure 


U(Y) 


u 


dx 


) 



(9) 


It is clear from equation (9) [and the figure] that the maximum velocity 
(in magnitude) occurs at the surface Y=l. Its value is given by U=l, 

The velocity is zero at Y=2/3 and reaches a maximum in the negative X 
direction at Y=l/3, the magnitude being |u| =1/3, 

This means that 


f = 

max 4p dx 


( 10 ) 


and that the actual x-directional fluid velocity will change by 4/3 this 
amount across the depth of the film. It is interesting to compare the 
order of magnitude of these velocities with the velocities expected for 
the bubbles migrating in the temperature gradient. A first approximation 
to the bubble migration velocity may be obtained from the results of 
Young, et al (1959) as 


^bubble 


4p dx 


( 11 ) 


Here, the variation of temperature with depth is ignored and , the 
interfacial tension between the liquid and the gas in the bubble, may 
be taken to be the same as the Interfacial tension between the liquid 
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I 


I 


I* 


I 


i 

f'ti 


i 


and the ambient medium. D is the bubble diameter. The negative sign indi- 
cates that the bubbles will migrate in the direction of the negative surface 
tension gradient (towards the hot end). Equations (10) and (11) imply 
serious consequences due to the surface-tension driven convection ^n the 
liquid film. It Is seen from these equations that the bubble velocities 
will, in general, be smaller than or on the same order of magnitude as 
the fluid velocity. Furthermore, the fluid velocity changes substan- 
tially from the top of the film to the bottom (these regions are 
indistinguishable to the optical system since all depths are in focus). 

Not knov/ing the vertical position of a given bubble, it will be virtually 
impossible to extract its migration velocity in the temperature field by 
subtracting out the fluid motion. It should also be pointed out that 
for a film depth of 0.25 mm and a bubble diameter of 0.2 mm, even if we 
knew the vertical position of the center of the bubble, it will be hard 
to simply subtract out the motion of the fluid since the effect of large 
variations of the fluid velocity over the bubble dimensions is complex. 
Intuitively, such subtraction might be acceptable if the bubble dimen- 
siua^ are less than about a tenth of the film thickness. Under these 
conditions, of course, the convection in the film induced by surface 
tension variation will play a dominating role on bubble migration. The 
only acceptable solution seems to be to eliminate such convection 
thereby eliminating the free surface. 

Example 

A set cf sample parameters are used below to calculate the 
expected order of magnitude of the velocities involved. Assume 

d ^ 0-025 cm 

2 2 
p ^ 10 poise (dynes sec/cm^) 

^ -0.06 dynes/cm 
^ -100 °C/cm 
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1 



f 




I 


Equation (10) gives the maximum velocity at the film surface as 

U -I ^ s- 3.75 X lO"^ cm/sec 
max 4|j dx 

For a bubble of diameter 0.005 cm, the migration velocity in the same 
temperature gradient, but in a stationary fluid, will be given by 
equation (11) as 


bubble^ 




12T 

4 m dx 


0.75 X 10 


cm/sec 
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I 
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I 


On the Stability of a i^olten Glass Film Heating Unifornily From Belov 

It is desirable to make at least some rough pre.liminaiy 
estimates of the crucial stability parameters in the molten gla^s film. 
This is done below. (The two kinds of instabilities to be anticipated 
are: (1) Surface-tension driven cellular convection, (2) Buoyancy- 

driven cellular convection (only on ground).) 

Marangoni Number 


Ma 


= - ip-) 

MT^ 


DAT 
U K 


( 1 ) 


AT = Temperature difference across the film depth 


dT 


Temp, gradient of surface tension 


- 0.06 


dynes 

cm 


D 

U 

K 




Film depth = 0.025 cm 

100 poise 

cm 

Thermal diffusivity % 0.01 cm^/sec 

of glass 


(rough estimate from 
Corning glass-properties 
booklet) 


The above set of parameters will yield a conservative 
estimate for the Marangoni number. The result is 


Ma ^ 0.0015 AT (2) 

where AT is in '‘C. So, even for AT on the order of lOO^'C 




Ma 0,15 


(3) 
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Now, for the onset of cellular convection driven by surface tension 
variations in the case of a J iquid film heated by an isothermal 
surface, learson gives 

Ma 'v 80 (A) 

c 

(It should be noted here that the present discussion is not related 

LJ the surface-tension driven flow examined in Appendix B which is 

caused by axial temperature gradients. These can only be avoided 

by placing a slide over the melt. The present calculations are 

concerned with Benard-type cellular convection in the absence of any 

temperature gradients along the heating surface.) Even though our 

Marangoni number is much smaller than Ma^ from Pearson's analysis, it is 

important to look farther since Pearson completely ignored gravitational 
18 

effects. Nield showed that in the presence of gravity, buoyancy 

effects reinforce surface tension effects — that is, instability occurs 

at smaller critical values of the Marangoni number than those predicted 

by Pearson. Nield's analytical results received experimental 

19 

verification a few years ago Palmer and Berg. Therefore, it is 
important to also investigate the Rayleigh number which characterizes 
buoyancy effects. 


Rayle ^ ,h Number (on ground) 


g 

D 

K 

V 


a 


Ra - 

K V 


acceleration due to gravity = 980 cm/sec 

film depth = 0.025 cm 

2 

thermal diffusivity of film ^ 0.01 cm /sec 

kinematic viscosity 
100 — 

cm* sec ^ 

= “ % 30 cm^/sec 

P ^ an_ 

^ 3 

cm 


thermal expansivity — 


increase in volume 
volume X 


10"^ (i/“c) 


Rough estimate 


(5) 
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The result is 


f 

» 

t 

N.I 


Ra 5 X 10 ^ AT (6) 

Again, with AT = 100®C, 

Ra ^ 5 X 10^^ (7) 

wbirh is extremely small compared to the critical Rayleigh Number 
Ra^ for the onset of buoyancy-driven convection calculated by Rayleigh* 

Ra^ % 1700 (8) 

The conclusion, therefore, is that neither buoyancy-driven 
instability nor Marangoni instability should cause concern in the 
glass film. Again, this conclusion refers only to cellular convection * 
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APPENDIX D 


Some Useful Information from Young et al 


8 


The final result in Young et al contains some typographical 
errors and should read; 

2 . . . 


37(2in~3n') ‘ ^ 


(9) 


for dimensional consistency and fat agreement with the Hadamard equation 
as However, whenwetried to derive this result starting from 

Eq. (5) to (8) in Young, ^ , we were unable to get exactly this 

result [or their Eq. (9)]. There appears to be a printing error 
in one or more of their starting equations. 

Assuming, temporarily, that the above equation is correct, 
it is quite reasonable to let g = 0 in it for simulating zero gravity. 
So, 

2uy'RT 


V = 
o 


3p(2y+3y ‘) 


in Young e_t £l . notation for this case. In view of some of the 
approximations made in arriving at this result, it is quite legitimate 
to make a few more y >> y* ; h >> h’ . This gives 


Here: 


V = -iRll 
2 M :tz 


V = velocity of bubble towards the hotter end 
R = bubble radius 

y = viscosity of the continuous phase fluid 


(10) 


lx - ^ ^ 


,)2 


dT 


2 


-e Temp gradient in the continuous phase 


Rate of change of surface 
tension with temperature. 
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